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classified and the effective potential is constructed in its most general form. 



transforms in the (3, 3) representation under the color SU(3) group. All 



the SU{3) invariant terms to the fourth-order of the order parameter are 



Ginzburg-Landau effective potential is studied for the order parameter that 




The conditions which stabilize the condensed phases are also obtained. As 



special cases and the Higgs phenomenon associated with symmetry breaking 



applications, the classification of the condensed phase is discussed in some 




is studied by introducing a coupling to the gauge boson. 
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I. INTRODUCTION 



Recently, quark-pair condensation phenomena (color-superconductor) have received 
much interest in the high-density nuclear physics [1,2]. 

In a high-density quark matter, quarks make a correlated Cooper pair due to the at- 
tractive one-gluon-exchange [3,6] or instanton-induced interactions [4]. A lot of studies have 
been done for physical properties of color superconductors on the basis of the QCD-inspired 
models [5-9]. 

From the symmetry point of view, the color-superconductor is characterized by the spon- 
taneous symmetry breaking of the color SUc(3) symmetry. Let's consider the quark field 
operator q a (r) with the color index a = 1,2,3 (other degrees of freedom, flavors and spins, 
are represented, e.g., as q = u^). They transform as a vector (the 3-representation) under 
the SU C (3) transformation. (Irreducible representations of 577(3) are denoted by their di- 
mensions, e.g. 1, 3, 3*, etc.) The color-superconductor is a long-range ordered state, the 
order parameter of which is a two-body one: ^ a ,p{f) = (q a (r)q'p(r)) . It is a second-rank 
tensor under SUc{3) with nine components, which makes a reducible 3 x 3-representation. 

Different from the systems with a scalar order parameter (metal superconductor, super- 
fluid 4 He, etc.), a variety of symmetry-breaking patterns shows up in the condensed states 
with the tensor order parameter. A typical example is the superfluid 3 He, where, due to the 
p-wave nature of the Cooper pair, the order parameter has a tensor-type spin structure, $jj 
= 1,2,3); it transforms as the 3 x 3-representation in the orbital and spin rotations 
G = 50^(3) <8> SOs(3) (also with the gauge U(l) symmetry). Corresponding to several sym- 
metry breaking patterns, many different condensed phases show up in the superfluid 3 He: 
A-phase (H = (5 2 x SO(3))/Z 2 ), B-phase (R B = 5i x 50 a (3)), and so on. As discussed 
later, such a phase structure should exist also for the S77(3)-3 x 3-type order parameter, 
and it is very interesting to study the structures for such a novel symmetry breaking case 
[7]. (It should be noted that the 3x3 representation of 577(3) is a complexification of that 
of 5*0(3), so that the color superconducting state can be considered as a kind of the group 
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extension for the superfluid 3 He.) 

The systems mentioned above are well described by the thermodynamic effective free 
energy that is constructed as a function of the order parameter up to the fourth-order, i.e., 
Ginzburg-Landau method [10]. This method incorporates the coherence nature of the pair 
condensation and describes the second order phase transition. It was also applied to the 
color-superconductor in some special cases [11]. 

In this paper, we discuss phase structures of the system with SU (3)-3 x 3 order parameter 
based on the Ginzburg-Landau method: the general form of the Ginzburg-Landau effective 
potential is presented to the fourth-order of the order parameter and analyzed in 

detail. As applications, the classification of the condensed phase and Higgs phenomenon are 
discussed in some special cases. 

II. STRUCTURE OF S77(3)-3 x 3 ORDER PARAMETER 

Let's consider the system with the symmetry group SU(3), where the order parameter 
is given by the 3-by-3 matrix : \I> = (^a,/?). For g G SU(3), it transforms as 

* — g&g, (i) 

where t g is a transposition of g. It makes a 3 x 3-representation of SU(3). 
The system should also be invariant under the U(l) gauge transformation: 

— ► e^W. (2) 

Because of the irreducible decomposition 3 x 3 = 3* + 6, the order parameter ^ = (ty a ,p) 
can be decomposed into the irreducible components: \I> = S + A, where S and A are the 
symmetric and antisymmetric components: 

S=^ + ^), A = ^(*-*tt). (3) 

Since the permutation symmetry for indices (a, (3) are conserved under the 577(3)- 
transformation, the S and A become basis for the subrepresentations of SU(3). By dimension 
counting, we find that the S and A correspond to the 6 and 3* representations. 
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From this consideration, symmetry breaking patterns for the condensed phases are clas- 
sified into four types: 



1) No-condensate A — S — 0, 

2) 3*-condensate A^O, S = 0, 

3) 6-condensate A = 0, S ^ 0, 
4) Mixed-condensate A^O, S ^ 0. 



(4) 
(5) 
(6) 
(7) 



Because of the 577(3) invariance of the system, for g e S77(3), order parameters (A, 5) 
and (gA t g, gS l g) represent the same condensed states, so that we can decompose the order 
parameters A and S ; 



A = e id GAoG, S = e ie GSjG, 



(8) 



where G G SU(3) and e ld e £7(1) are the Goldstone degrees of freedom which can be 
eliminated by the 577(3) and £7(1) gauge transformations. Remaining parts 7l and 5 in 
eq. (8) determine essentially scales of the symmetry breaking. As proved in appendix A, 
they can be parametrized as 

/ 03 -02 \ M \ 
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(9) 



where di,2,3 are real and 01,2,3 are complex parameters: 0, = |0j|e* e \ We take the 
parametrization in eqs. (8) and (9) as a standard form and call di and 0« normal coor- 
dinates of the 3x3 order parameter. 



III. EFFECTIVE POTENTIAL FOR S77(3)-3 x 3 ORDER PARAMETER 

In the Ginzburg-Landau theory, the effective potential V e s for the order parameter \1/ 
plays a central role. It should be a invariant function of \I> for the SU C (3) and U(l) gauge 
transformations, (1) and (2). For a minimal Ginzburg-Landau model, it is necessary to 
construct V e & to the fourth-order of 
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The lowest order contributions to V e s are the second-order terms in \1> and there exist 
such two terms constructed from A and S: 

V$ = ai TiA*A + a 2 TiS*S. (10) 

All other second-order terms can be shown to vanish or to be transformed into the above 
form (e.g., TrA*S = due to the antisymmetry of A*S.) 

In a similar way, we can construct a general fourth-order potential v}g\ 

V$ = b 1 (TrA*A) 2 + b 2 (TrS*S) 2 + b 3 TrS*SS*S 

+ b 4 TrA*ATrS*S + b 5 TrA*AS*S + b+(TrA*SA*S + c.c) + ib-(TrA*SA*S - c.c). (11) 

In the same way as for V^ 2 \ any other fourth-order terms can be transformed into the above 
form: e.g., (TrAM) 2 = 2TrA*AA*A. 

Using eq. (8), the potentials V^' 4 \ eqs. (10) and (11), can be rewritten with the normal 
coordinates (di,\4>i\,9i) defined in eq. (9). For simplification, we introduce the variables 
R = ^J\(f)i\ 2 + |02| 2 + |03| 2 j D = \Jd\ + d\ + d\ and the normalized coordinates: r; = \4>i\/R 
and Ci = di/D, then eqs. (10) and (11) become 

K ( ff 2+4) = 9rR 2 + 9rR A + 9dD 2 + (g' D - h D x)D A + {g" D + h" D u)R 2 D 2 1 (12) 
where the coefficients gi are defined by linear combinations of a, and b{ in eqs. (10) and (11): 

g R = -2a u g' R = 4b 1 , g D = a 2 , g' D = b 2 + b 3 , 

h D = 2b 3 , & = -2b 4 -b 5 , h" D = b b . (13) 

The variables x and u in eq. (12) are defined by 

x = e\e\ + e 2 e 2 + e 2 e 2 , u = j x r\ + f 2 r 2 2 + / 3 r 2 , (14) 

where 

/i = e 2 + csin(2^i + 5)e 2 e 3 , f 2 = e\ + csin(2# 2 + 5)e 3 e X) f 3 = e 2 + csin(2^ 3 + 5)eie 2 . 

(15) 

The constants c and 5 in eq. (15) are obtained from the coefficients b± in eq. (11) as c = 
-4^/6+ + &V&5 and tan5 = 
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IV. DOMAIN OF VARIABLES 



Now we obtain the general form of the effective potential to the fourth-order, which is a 
function of the variables (R,D,x,u). Let's consider the domain of these parameters. 

First, the variables R and D are defined as the lengths of |0j| and di, so that R > and 
D > 0. 

We introduce new variables: 

0<pi = e?<l, 0<qi = rj <1, -1 < \ = sin(^ + 8) < 1, (i = 1, 2, 3) (16) 

which satisfies the constraints: P1+P2+P3 = 1 and q\ + q 2 + q 3 = 1. Geometrically, possible 
ranges in the (pi,P2,P3) or (gi, 52,93) space become the two-dimensional region inside the 
equilateral triangle with the vertices (1,0,0), (0, 1,0) and (0,0, 1). Using these variables, 
the variables x and u are represented by 

x = P2P3 + PsPi + P2P3 = ^ {l - (pl +pl +pI)} , u = /igi + f 2 q 2 + fsqs, (17) 

where f 1=Pl + c\ lv /p 2 p 3 . 

Let's consider the domain of x — pip 2 + P2P3 + p 3 Pi- It is clear that x > because 
Pi,2,3 > 0. Using the inequalities (pi + Pj) 2 > 4=PiPj, we obtain 

1 , , , x 2 (Pl+P 2 ) 2 + (P2+P 3 ) 2 + (P3+Pl) 2 , . 

l = ( Pl +p 2 +p 3 ) = ha; > — + a; = 3a;. (is) 

Thus, the variable x takes the values < x < 1/3. 

The domain of tt = /i^i + /2<?2 + /3<?3 depends on the value of x. Here, we summarize 
the derivation: 

1. Minimum /maximum of u for fixed f\. Because < qi < 1 with q x + q 2 + q 3 = 1, the 
domain of u becomes min{/i, / 2 , h} < u < max{/i, f 2 , f 3 }. 

2. Minimum/maximum of f\ — p\ + cAi^/p^pi /or /ised p^. Without loss of generality, we 
can assume c > because — 1 < Ai < 1. For fixed pi, the variable f\ takes a minimum 
value pi — Cy/p2P3 at Ai = —1 and a maximum value pi + c^/p^pi at Ai = 1. Because 
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fi are defined symmetrically for the index (i = 1, 2, 3), similar results are obtained for 
f 2 and / 3 . 

3. Range of u for fixed pi. Combining the above results 1. and 2. , we obtain the possible 
range of u for fixed p{. /{" < u < /{*", where 

fr =Pi- c^/p~m, ft = pi + cjpipi. (19) 

4. Finally the domain of u for the fixed value of x is obtained by min {/f, f 2 ~, /r} ^ M ^ 
max /2", /jj 1 "} where min and max are for < pi < 1 with pi +p 2 +P3 — 1. 

Thus, the evaluation of the domain of u has been reduced to those of min(/_) and max(/ + ), 
which are somewhat complicated because they depend not only on the coefficient c but also 
on the variable x. The calculations for them are given in appendix B, and we show only the 
final results for the domain of (x, u) below. 

The results are given for (0 < c < 2) and (2 < c), which are respectively shown in 
figure 1(a) and (b). These figures show upper and lower bounds of u for fixed x. 

a : < c < 2 

2 — c 1 + c 

Upper bound of u : yl — 3a; H (20) 

3 3 



Lower bound of u : < 



-2±VT=3£+i=£ (*<*<§) 



The segment AB in figure 1(a) corresponds to eq. (20), and OC and CD to two equations 
in (21). 

b : 2 < c 



Upper bound of w : < 



-c^ (0 < re < i) 

Lower bound of u : ^ (23) 



The segments AE,EF and FB in figure 1(b) correspond to three divisions in eq. (22), and 
OC and CD to two equations in (23). 



V. STABILITY CONDITION 

For the system with the effective potential V e g to be physically meaningful, it should be 
positive definite when the variables R or D take infinitely large value. (If not, any states 
with finite R and D become unstable.) 

This stability (positive definite) condition gives constraints on the parameters and hi 
in eq. (12) and we derive them in this section. First, we consider special cases: 

1. When D = and R — > oo, eq. (12) become V e g ~ fifcR 4 

2. When R = x = and D — > oo, we obtain V c g ~ g^-D 4 . 

From the positive definineness of V^r in these limits, the first two constraints can be obtained: 

g' D >0, g' R >0. (24) 

Hereafter, we consider the cases of g' D > and g' R > only. (We can treat the cases that 
g' D = or g' R = with taking the limit of those in g' D > and g' R > 0.) 

We replace the variables i? and by the scaled ones: R — > i?/ (g 1 ^)^ and D ^ D / (g' D )^, 
then eq. (12) becomes 

Kfr = -Gi? 2 + i? 4 - FD 2 + (1 + F'a;)^ 4 + 2(// + H'u)R 2 D 2 , (25) 

where scaled parameters are defined by 

G = _J^ F F' = -^, H=^^, # = (26) 

V^k V^d fe 2 y9D9R 2 \J9'd9'r 

For large values of i? 2 and -D 2 , eq. (25) behaves as a linear quadratic form of them for 
any fixed values of (x,u): V eS ~ R 4 + (1 + F'x)D 4 + 2(if + H'u)R 2 D 2 . The Sylvester's 
theorem should give the positive definite condition for it [12]: 
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a) 1 + F'x-(H + H'u) 2 > 

or 

b) 1 + F'x > and H + H'u > 0, (27) 

for any values of (x,u), the domain of which have been shown in figure 1(a) and (b). The 
second condition b) in eq. (27) should be required because R 2 , D 2 > 0. 

The condition (27) gives constraints for the possible values of parameters F', H and H', 
which can be obtained analytically. We sketch the derivation of them in Appendix C, and 
only show the results here. 

Stability conditions in (H, F') space for fixed c and H' are classified at first into two 
types for (0 < c < 2) and (2 < c), and then into four types by the values of H' in each of 
the first two (which correspond to four subscripts I~IV in figure 2 and figure 3), which are 
shown in figure 2(a)~(d) for < c < 2, in figure 3(a)~(d) for 2 < c. In the following lists, 
we obtain complete stability conditions for the coefficients F', H, H' and c in eq. (25). 

A. < c < 2 case 

a : H' < -^r c (figure 2(a)) 

F' > -3, (-H' - 1 < H) (28) 

This critical point {—H' — 1, —3) corresponds to the point A in figure 2(a). 
b : -^r c < H' < (figure 2(b)) 

F' > k(c), (-H' -1<H< -^-^ H ') ~AB (29) 
F' > -3, (-^-^H' < H^j B~^o~ (30) 



where k(c) 



II + H') (H + ^H>) -1-J{(H + H'f - 1} { (H + ^H'f - l} 



c: < H' < 4^ (figure 2(c)) 



F' > ^Qt, (-1 < H < ^o~A 

F' > (2H - cH') 2 - 4, {H X <H< F 2 (c)) AB 
F' > fc(-c), (F 2 (c) <H< ^#') BC 



2 



cH' - 1 



F' > -3, \-^-H'<Hj Coo 

where H, ee F 2 (c) = (^W^) 2 ^ . 

^<F' (figure 2(d)) 

F ' ^ (-1 < F < #i) ^4 

F' > (2F - cF') 2 - 4, [Hi < H < CH '~ 1 ) AB 



F' > -3, < if Foo 



B. 2 < c case 



a: H'<—^ (figure 3(a)) 

F' > (2F + cH'f - 4, - 1 < F < 



l + cH r 



1 + cH' 



F > k(c), ( F 2 (-c) < F < —H' ) FC 

1 + C rr/ 



AF 



F' > -3, y < Hj Boo 

b: < F' < (figure 3(b)) 



F' >{2H + cH'f -A, (-H' - 1 < H < H 2 (-c)) AB 

1 + C rw 



F' > -3, ( — F' < F ) Coo 

c: < H' < 4^ (figure 3(c)) 
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F' > 



(cH 



/\2 



1 — if 2 ' 
F' > (2H - cH') 2 - 4, 

F' > fc(-c), 
F' > -3, 



d: ^ < H' (figure 3(d)) 



F> > {CH ' )2 



1-H 2 ' 
F' > (2H - cH') 2 - 4, 

F' > -3, 



(-1 < H < 

(Hi < H < H 2 (c)) 



H 2 (c) < H < - — -H' 



c- 1 



H' < H 



(-1<H< 

cH' -V 



Hi<H < 



'cH' - 1 



< H 



ooA 



AB 



BC 



Coo 



ooA 



AB 



Boo 



(43) 
(44) 
(45) 
(46) 



(47) 
(48) 

(49) 



VI. APPLICATIONS 

In this section we show some applications to our model. As a first application, we obtain 
the phase diagrams for the two special sets of the coefficients (26). As a second application, 
we construct a kinetic term by introducing a coupling to the gauge boson and study the 
Higgs phenomenon in the three characteristic phases. 



A. Phase diagrams for some cases 

1. The angular parameter (x,u) independent potential 

We study the condesed phases in the case of (x, w)-independent potential, F' — and 
H' = in eq. (25). The effective potential becomes a function of R and D with three 
parameters (G,F,H): 

V(R, D) = -GR 2 + R 4 - FD 2 + D A + 2HR 2 D 2 . (50) 

The stability conditions are given as in the previous section: 
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1 - H 2 > or H > 0. (51) 

Four kinds of condensesd phase (as extremum points of the effective potential) are ob- 
tained: 1) R = D = , 2) #0 ^ 0, D = , 3) R = 0, D ^ , 4) R ^ 0, 
Do 7^ corresponding to eqs. (4)- (7), and the parameter regions for each phase are shown 
in Table 1 with the conditions (51). Referring to Table 1, the phase diagram on (F, G) 
space for fixed H is obtained in figure 4 by comparing V(Ro, D ) in the ranges where more 
than two kinds of extrema exist. At H——1, the mixed phase emerges above the line AB 
(G——F) in figure 4(a), and as H increases to 1 the two branches of this line merges into 
the line OA (G = F) in figure 4(c). When H exceeds 1 the mixed phase vanishes. 

2. The potential for the symmetric phase 

Next we study the case in which only the symmetric 6 phase emerges, i.e., 
9r = 9'r = 9d = ^d = ^ m ec l- (12). In this case the effective potentail is given by 

V(x, D) = g D D 2 + (g' D - h D x)D A . (52) 

From < x < |, we obtained the stability condition by 

h D <3g' D i 9d> q - ( 53 ) 

This condition coincides with the previous results in F>— 3 and H'—H— 0, e.g., in eqs. (29) 
and (30). From eq. (52) the four kinds of condensed phase are obtained: 1) D =0 , x—0, 2) 
-DoT^O, x—0, 3) -DoT^O, 4) D ^0, The parameter regions and the form of the order 

parameter So in eq. (9) for each phase are given in Table 2. The phase diagram in (F, G) 
space for g_D>0 is shown in figure 5. 

It is noticed that the linear dependence of the effective potential (52) on x is the reason 
why the phase characterized by, e.g., So—(di, d,2, 0) does not emerge. 
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B. Higgs phenomenon 

Higgs phenomenon is the one where gauge bosons acquire finite effective masses in the 
broken-symmetry phase. In the minimal coupling scheme, the gauge boson fields are 
introduced through the kinetic terms of the order parameter. To the second-order, it is 
given by 

K kmeUc = -k A Ti(D^A)*D^A + ksTr(D IJ/ S)*D^S, (54) 

where the coefficients k A and ks are positive constants. The covariant derivative in 
eq. (54) is defined by 

= dfl - ig^A^ - ig^A;, (55) 

where g and A a= i^ 8 are the gauge coupling constant and the Gell-Mann matrices [14]. 

The effective masses of the gauge bosons are obtained from the second-order terms of 
the gauge boson fields in the kinetic terms (54). 



1. The 3* condensed phase 



In the case of the 3* condensed phase we consider an order parameter 

/0 -a 0\ 

^o= a , (56) 
v 0/ 

where a is a real constant. Note that this form is general enough since other forms of the 
3* phase is obtained by a simple transformation to ^ as in eq. (8). In this phase SU C (3) 
symmetry breaks into SU(2). 

By substituting \I> into eq. (54) the mass terms are obtained: m 2 v Tiv^v and m?Tr0 2 



where m v =\/2kA\a\g, m^=y ^k A \a\g, v=^(X 4 + i\ 5 )(A* + iA^) +§(A 6 + i\ 7 )(A^ + iAjJ and 
(p=\ 8 A 8 ^. The gauge fields v and <f> are different multiplets of the residual symmetry SU(2). 
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This phase is the one which is featured in most microscopic models. Using the values of 
the gauge coupling constant g and the gap \a\ caluculated in models [4], the effective masses 
are estimated: m v =J^ms ~ 100 MeV. 



2. The 6 condensed phase 

We select an order parameter proportional to the identity matrix which is discussed in 
the previous application, see Table 2 : 

fa 0\ 



(57) 



a 
\0 a) 

In this phase SU (3) symmetry breaks into SO(3). Since the gauge boson field A^(=\ a A a ^) is 
transformed as direct products of two triplets under the SO (3) group, they are decomposed 
into the direct sum of triplet (anti-symmetric part) and quintet (symmetric part): 



A„ = Va + T B , (V a = £ A a A«), (T 5 = £ \ a Al). 

a=2,5,7 o=l,3,4,6,8 



(58) 



Since T 5 is assigned to the generators of the broken symmetry SU{3)/ SO(3), the mass term 
is accompanied with the quintet term: m^Tr(T 5 ) 2 , where m T =^/Aks\a\g. 

3. The Mixed condensed phase 

We study a special case for the mixed phase characterized by the order parameter: 

/0 0\ 

= . (59) 
\a 0/ 

In this phase SU(3) symmetry breaks into U(l) which is generated by H ^(t) =e *( A 3+v / 3A s )t _ 
Since in this case the multiplets of the gauge field are not trivial, we investigate a dependence 
of the gauge field on the infinitesimal transformation of H(t). 
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( 1 2\ 
-1 1 

v-2 -l oy 



(60) 



Noticing that (A 3 + y/SXs) and (A 3 — ^A 8 ) are orthogonal each other, the gauge field are 
decomposed into the multiplets: 

/0 f v\ 



(A 3 - ^ A sV + (A 3 + >/3A 8 y , 



£ = A}~n ~ iA 2 ^, r] = A^ — iA 5 ^, ( = A 6 ^ — iA 7 ^ 
<p ee (3Al - V3AD/4, tf/ ee (Al + V3AD/4. 



The mass terms are summerized as follows: 



(61) 



m|re = ^(k A + k s )\a\ 2 g 2 C^, rn^rj = 2k s \a\ 2 g 2 r ] *r ] 

1 16 

mlCC = ^ + k s )\a\ 2 g 2 CC, m^ 2 = —(k A + k s )\a\ 2 g 2 V 2 



(62) 



In the same way more complicated cases can also be studied. 



VII. SUMMARY 

We have derived the general form of the Co\oi-SU (3)-Ginzburg-Landau Effective Poten- 
tial for the order parameter with 3x3 symmetry to the fourth-order and determined the 
range of coefficients which stabilizes the vacuum of the system. Corresponding to the irre- 
ducible components of the order parameter 3x3=3*+6, we found four kinds of phases: non 
condensed phase, anti-symmetric (3*) phase, symmetric (6) phase and mixed (3* +6) phase. 

As applications, we obtained the phase diagrams in the two special cases and studied 
the Higgs phenomenon in some characteristic phases by gauging this model. 

The full classification of the broken vacuum are analytically possible including the angular 
parameter x and u, and will be discussed in different paper. 
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APPENDIX A: NORMAL COORDINATE DECOMPOSITION 



In this appendix, we show the decomposition of the 3-by-3 symmetric matrix S — f S in 
eq. (9). The proof is a variation of that used in the group representation theory [13], the 
SU (3)-chiral nonlinear sigma model [14] and the Kobayashi-Maskawa theory of quark mass 
matrix [15]. 

We take the subsidery hermite matrix H = S*S = S^S. Its eigenvalue Aj and eigenvectors 
are obtained from the eigenequation: ifej = A^e, (i = 1, 2, 3). 

In this appendix, we consider only the case of no degeneracy: Aj ^ Xj (i ^ j). (The 
extension to degenerate cases should be made by analytic continuation. ) 

Because of the hermiticity of H, the eigenvectors e, can be taken to be orthogonal unit 
vectors (^e^ej = <^), and the eigenvalues A« take real positive values, (because of H = S^S). 
In matrix form, the eigenequation becomes 

/Ai \ 



U ] HU = L = 



A 2 

Vo o a 3 ; 



(Al) 



where U = (ei,e 2 ,e 3 ) E SU(3). 
The eigenvector ej satisfies 



H*Se t = SS*Sei = SH ei = \ t Se t . 



(A2) 



so that we obtain Sei = «je* (because H*e* = A«e*). ai = \ai\e l9i are complex numbers. In 
matrix form, it becomes 

( ai \ 



*usu = z = 



a 2 
V a 3 J 



The polar decomposition of Z becomes 

/ Ictil \ / e idl \ 



A = Re 



iT 



o H o 

V \a 3 \J 



e i(>2 
V e i03 J 
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, T 



(Q\ 0\ 

o e 2 o 

VO 9 3 J 



(A3) 



(A4) 



Because of A^e, = Hei = S*Sei = ctiS*e* = cti(Sei) = |aj| 2 ej, we obtain A« = \a\i and 
also R = \fh. The polar decomposition of A becomes A = \fZe lT . 

Because the angular part T is diagonal, it can be expanded by the 3-by-3 unit matrix 
I and the diagonal Gell-Mann matrices A 3j8 : T = 91 + 3 A 3 + 4>s\s- If we define K = 
efOMa+^As^ then we obtain e iT = Ke i6 IK. Using the first equation in (A4), the matrix S 
becomes 

S = U*ZU = U*^e iT U^ = e w U*KVLKU ] = e^G^G, (A5) 

where G = U*K and f G = t KU jf = KW . This is just the decomposition of S in eq. (8) with 
S = VI. 

If we define the matrix A by A = GAqG, then A also becomes antisymmetric and is 
generally represented as the first equation in eq. (9). 

APPENDIX B: MAXIMUM AND MINIMUM OF / 

Instead of pi with p± +p 2 +P3 — 1, we take independent variables (s, t): 

S = ^=(~Pl+P2), ^-^(Pl+Pl-I)- ( B1 ) 

The possible range of (s, t) is determined from the condition: < Pi < 1 and P1+P2+P3 = 1- 



Using (Bl), it is found to be the inside of the equilateral triangle with vertices (0, y2/3) 
and (±l/\/2, — 1/a/6) in the (s,t) plane. Since the function f(s,t) is an even function of s, 
f(s,t) = f(—s,t), we may consider a half (s > 0) of the triangular region (figure 6) to find 
its maximum/minimum. 

For a fixed value of x, a trajectory of x — P2P3+P3P1+P1P2 = — (s 2 + t 2 )/2 + l/3 in (s,t)- 



plane become a circle with the center (0,0) and the radius yj (2/3) — 2x (figure 6). From 
the aspect of intersections between the circle and the triangular region, two cases should be 
discriminated: 
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1. When < x < 1/4, the circle intersects with three sides of the triangle (figure 6(a)), 
and the two separate arcs (B Q B 1 and B 2 B 3 in figure 6(a)) are included in the half 
triangular region. 

2. When 1/4 < x < 1/3, the whole circle is included in the triangular region (figure 6(b)). 
The intersections with the s-axis (t = 0) are denoted by B and _B 4 . 

The (s, t)-coordinates of the above four points -Bq~4 are given by 



B,= [0,-^-2x\. (B4) 
Let's consider the maximum value of f£ for fixed values of c and x. Using s = 



V6 - 9t 2 - lSx/3 and (B4), / 3 + in eq. (19) becomes 

~ (B5) 



f+(t) = ^ V6t + 1 + c\/6t 2 - ^ + 9^-2 
where we take t as an independent variable and (x, c) as any fixed parameters. Solving 



the equation [f^it)}' = 0, we obtain one local extremum at B M = (yo — 9t 2 M — 18x/3, tj^) 
where 

1 



tu — 



2VQ 



(B6) 



1 -c 2 

for some values of (x, c). The parameter range where the extremum of f£ exists (at Bm) can 
be read off from the condition that t M in (B6) should take a real value: i.e., (Ax — l)(c 2 — 1) > 
0. Solving this inequality, we obtain two separated regions: a)c< 1 & x < 1/4 and b) c > 1 
& x > 1/4 (Shaded regions in figure 7) where the function f% has an extremum. Judging 
from the signature of [/^"(iw"), we can find that /^(^m) is a local maximum in region a) 
and minimum in region b). For tM have to be in the permissible region shown in figures 6(a) 
and (b), we obtain further condition: 
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1. When x < 1/4 and c < 1, t M have to be on the arc B 2 B 3 (figure 6(a)). Thus, we 
obtain x > c 2 /4 (OD in figure 7). 

2. When x > 1/4, t M have to be on the half circle B B 4 (figure 6(b)). Thus, we obtain 
c < 2 or x < c(5c - 4)/[4(2c - l) 2 ] (AS in figure 7), 

Summarizing the above discussions about the aspects of extrema of f 3 , we can classify the 
five regions in the (c, x)-plane (figure 7). 

Let's analyze maximum points in each regions. Candidates of the absolute maximum 
point of /| are the end points -Bq~4 where 



Comparing f+(Bi) in eqs. (B7) and (B8), we can find maximum points in each region in 
figure 7: 

1. The maximum is at B in regions 1, 2, apart of 3 (c < 2 or x < (2c— 1) 2 /[4(1 — c+c 2 ) 2 ] 
EF in figure 8 ), and a part of 4 (c < 2). 

2. The maximum is at B 3 in a remaining part of region 3. 

3. The maximum is at B 4 in region 5 and a remaining part of 4. 

The results are summarized in figure 8. 

In a similar way, the minimum value of fs(t) can be calculated: 



fl(B ) = i±£ + ^yr^, 4( J Bi) = - + - v / r^, 



fl(B 2 ) = \ - ±VT=Ix~, fl(B 3 ) = cVx~, 



(B7) 



or the local maximum at M (if exist): 




(B8) 




(B9) 
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APPENDIX C: DERIVATION OF STABILITY CONDITIONS 



In this appendix, we sketche the derivation of the stability condition: a parameter range 
of (F',H,H') that satisfies (27). 

Instead of (x,u), we take new variables (X, U) defined by 

X = l + F'x, U = H + H'u. (CI) 

Using them, the conditions in (27) become simple: 

a) X > U 2 , or b) X > and U > 0. (C2) 

In the (X, C/)-plane, the range that satisfies these conditions is represented by the hatched 
area in figure 9. We denote it by the region E. 

The possible range of variables (X, U), which is translated from the one of (x,u) (given 
in figure 1), depends on parameters (F', H', c) and we denote it by A(F', H', c). In figure 9, 
as an example, we show the A(F', H', c) as the cross-hatched region for some values of 
(F', H', c) that satisfies F' > 0, H' > and c > 2. In (X, C/)-plane, the stability condition 
(27) become equivalent with the geometrical condition: 

A(F', H', c) C E. (C3) 

The boundary of A(F', H', c) and E are patchworks of line or parabolic segments, so that 
the condition (C3) can be solved algebraically. 

As an example, let's consider the case that c < 2, F' > and H' > (figure 9). In this 
case, eq. (C3) can be attributed into whether the boundary segment ooGO for E intersects 
that of A, AB, or not. Both segments, ooGO and AB, have parabolic shapes represented 
by X = U 2 (ooGO) and X = [F' / (cH') 2 ](U — H) 2 + 1 (AB). It can be judged from the 
condition that the edge points, A = (1,H) and B = + 1, H — are included in E; it 
gives two necessary conditions for parameters: 

-1<H, (H- C -pj 2 <^ + l. (C4) 
20 



In case that F'/(cH') 2 < 1, the slope of ooGO is more gradual than that of AB, and eqs. 
(C4) are sufficient condition. However, in the case that F'/(cH') 2 > 1, an extra condition 
should be added in order that two parabolic segments have no intersections: 

(cH'f + F\H 2 - 1) < 0. (C5) 

Eqs. (C4) and (C5) give the necessary and sufficient stability condition in that case: 

(cH') 2 

F'>j—^, F' > (2H - cH') 2 - 4, -1<H, (C6) 

which is shown as the hatched region in figure 2. 
Other cases can be treated in a similar way. 
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figure 1. Range of variables x and u. a) for < c < 2, b) for 2 < c. Details are given in 
the text. 

figure 2. Stability region in (H,F') space for (0 < c < 2) (the shaded area), figure 2(a) is 
for H' < -£- c , (b) for < H' < 0, (c) for < H' < ^ and (d) for ^ < H' . 

figure 3. Stability region in (H,F r ) space for (2 < c). figure 3(a) is for H' < —-^i (b) 
for < H' < 0, (c) for < < ^ and (d) for ^ < H' . 

figure 4. (a) The phase diagram for the effective potential (50) in (F, G) space when H——1. 
AB : G=—F. The regions with vertical stripes, with horizontal stripes and with cross 
hatch show 3*, 6 and the Mixed condensed phases. The remaining part (G < and 
F < 0) shows the Non condensed phase, (b) The same, but H=l/2. OA : G=2F. 
OB : G=F/2. (c) The same, but H=l. OA : G=F. 

figure 5. The phase diagram for the effective potential (52) in (g' D , ho) space when #d<0. 
AB : h D =3g' D . The regions I and II corresponds to the phases 2) and 4) in Table 2, 
but the region just on the line /id=0 (^>0) corresponds to the phase 3). 

figure 6. Range of variables s and t is shown as inside of the half triangle ACD. Trajectories 
of (s, t) for fixed x are given for < x < \ by two separate arcs B$Bi and .B2-B3 in 
(a) and for \ < x < | by a half circle B B 4 in (b). 

figure 7. The five regions in the (c, x)-plane according to the aspects of tu- tu is defined 
in the shaded regions. In addition, in the region 1 and 4 tu exists on the arc B 2 B 3 
(figure 6(a)), and on the half circle -B0-B4 (figure 6(b)). 

figure 8. The three regions in the (c, x)-plane separated by the maximum values of /+, 
which are f+(B ) , f+(B 3 ) and f\{B A ) in the region I, II and III. 

figure 9. An illustration for the stability condition A(F', H', c) C S in the case c < 2, F 1 
> and H' > 0. 
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Table 1 : Condensed phases and their parameter ranges. 





i?o 


D 2 


V(R , D ) 


the ranges where Ro, D are defined 


Non cond. 













G<0, F<0 


3* cond. 


G/2 





-G 2 /4 




G>0 


6 cond. 





F/2 


-F 2 /4 




F > 

G > HF and G > F/H, if -1 < H < 


Mix cond. 


G-HF 
2(1-H 2 ) 


F—HG 
2{1-H 2 ) 


-G 2 -F 2 +2FGH 
4(1-H 2 ) 


< 


G > HF and G < F/H, if < H < 1 
G < i/F and G > F/H, if 1 < H 
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Table 2 : Forms of the order parameter and their parameter ranges 



(Dl x) 


So = (di,d 2 ,d 3 ) 


the parameter region 


1) (0, 0) 

2) (-ft- 0) 

3) (-ft, x) 

4 ) ( 2( s ' n -ft D /3)' 3) 


di = d 2 = d 3 = 
di 0, dtfi = 
d? + d 2 2 + d§ = D 2 
d\ = d 2 = d 3 ^ 


9d>0 

g D <0, h D <0 
g D < 0, h D = 
#d < 0, < /i D < 3^ 
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